mm m 

INTEGRATIVE NEUROSCIENCE 



ORIGINAL RESEARCH ARTICLE 

published: 13 March 2013 
doi: 10.3389/fnint.2013.00013 




Axonal and glial microstructural information obtained 
with diffusion-weighted magnetic resonance 
spectroscopy at 7T 

Itamar Ronen*, Ece Ercan and Andrew Webb 

C.J. Gorter Center for High Field MRI, Department of Radiology, Leiden University Medical Center, Leiden, Netherlands 



Edited by: 

Alpay Ozcan, Virginia Polytechnic 
and State University, USA 

Reviewed by: 

Christopher D. Kroenke, Oregon 
Health and Science University, USA 
Dost Ongur, McLean Hospital, USA 
Julien Valette, Commissariat a 
I'Energie Atomique, France 

Correspondence: 

Itamar Ronen, Department of 
Radiology, Leiden University 
Medical Center, Albinusdreef 2, 
Leiden, 2333 ZA, Netherlands, 
e-mail: i.ronen@lumc.nl 



Diffusion-weighted magnetic resonance spectroscopy (DWS) offers unique access to 
compartment-specific microstructural information on tissue, and potentially sensitive 
detection of compartment-specific changes in disease. The specificity of DWS is, however, 
offset by its relative low sensitivity, intrinsic to all MRS-based methods, and further 
exacerbated by the signal loss due to the diffusion weighting and long echo times. In 
this work we first provide an experimental example for the type of compartment-specific 
information that can be obtained with DWS from a small volume of interest (VOI) in brain 
white matter. We then propose and discuss a strategy for the analysis of DWS data, which 
includes the use of models of diffusion in compartments with simple geometries. We 
conclude with a broader discussion of the potential role of DWS in the characterization 
of tissue microstructure and the complementarity of DWS with less-specific but more 
sensitive microstructural tools such as diffusion tensor imaging. 

Keywords: tissue microstructure, magnetic resonance spectroscopy, high field MRI, diffusion-weighted imaging, 
white matter, N-acetylaspartate, choline, creatine 



INTRODUCTION 

The characteristic dimensions of microstructural features of liv- 
ing tissue are several orders of magnitude smaller than the typ- 
ical spatial resolution obtained by magnetic resonance imaging 
(MRI). Yet, MRI has become a ubiquitous tool for investigating 
tissue microstructure in vivo, thanks to the direct relationship 
between the motional properties of the nuclei that generate the 
MR signal and the properties of the NMR signal that ensues. In 
particular, self-diffusion of spins in a magnetic field gradient gen- 
erates a signal loss that is controlled by the characteristics of the 
gradients as well as by geometric factors that hinder or restrict 
the diffusion process over time. Incorporation of pulsed magnetic 
field gradients in MR pulse sequences led to the most robust man- 
ner of sensitizing the MR signal to self-diffusion (Stejskal and 
Tanner, 1965) and has since been frequently used in the inves- 
tigation of numerous aspects of tissue microstructure, including 
characterization of disease processes (Le Bihan et al., 2001), 
exploration of white matter connectivity (Conturo et al, 1999; 
Melhem et al., 2002), and quantitative evaluation of microstruc- 
tural features such as axonal diameter and myelination (Nomura 
et al, 1994; Neil et al., 1998; Assaf et al, 2008). 

The explanatory power of diffusion-weighted (DW) MRI 
methods is limited, however, by the lack of specificity of the sig- 
nal source, namely, the hydrogen protons in water molecules. The 
ubiquity of water in living tissue (~65% in volume) accounts 
for the high sensitivity of MRI when compared to MR of other 
tissue constituents, but the lack of compartmental-specificity 
complicates the interpretation of DW-MRI results. For example, 
competing disease processes that simultaneously act on several 



tissue compartments may have similar, or opposed effects on dif- 
fusion properties of water in those compartments. This, in turn, 
generates a high degree of ambiguity in the interpretation of 
the DW-MRI results, or in other instances dilutes the informa- 
tion related to the specific process of interest. Examples include 
DW-MRI of stroke, where the decrease in apparent diffusion coef- 
ficient (ADC) in the ischemic region can be attributed to several 
competing effects that influence the ADC: increase in restric- 
tion in the diffusion of water in the extracellular compartment, 
influx of water into the restricted intracellular compartment, and 
the cessation of active transport (cytosolic flux) in the affected 
cells. In the case of multiple sclerosis (MS) and other neuro- 
logical disorders, damage to myelin may increase water ADC 
and reduce the fractional anisotropy (FA) of the diffusion ten- 
sor, but neuronal/axonal damage may reduce overall ADC. It is 
thus desirable to devise a MR-based modality with an intrin- 
sic ability to provide compartment-specific information about 
tissue microstructure. Several approaches to achieve this goal 
have been suggested, including high-level modeling of DW-MRI 
acquired at a variety of diffusion-weighting values and direc- 
tions (Assaf and Basser, 2005; Zhang et al., 2012). A more direct 
way to tackle the issue of compartment-specificity is to measure 
the diffusion properties of endogenous spin species that reside 
in specific tissue compartments and are concentrated enough to 
yield a visible MR signal with reasonable signal-to-noise ratio 
(SNR). In neural tissue, these are intracellular metabolites in 
the mM concentration range and give rise to distinct spectro- 
scopic patterns in a magnetic resonance spectroscopy (MRS) 
experiment. Some metabolites are almost entirely cell-specific, 
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e.g., N-acetylaspartate (NAA) and the excitatory neurotrans- 
mitter glutamate (Glu), which reside both almost exclusively 
in neurons/axons, and myo-inositol (ml) and glutamine (Gin), 
which are thought to primarily reside in astrocytes (Brand et al., 
1993; Choi et al, 2007). Other intracellular metabolites, e.g., cre- 
atine and phosphocreatine (Cr and PCr) are less-specific and 
have a well-defined role in cellular metabolism, while others are 
preferential to a specific cell type but not exclusively so, e.g., 
choline-based molecules, which are believed to be more con- 
centrated in glia than in neurons at a ratio of approximately 
3:1 (Choi et al, 2007). Given the low concentration of these 
metabolites, the sensitivity of their detection is much lower 
than that of water, and thus the spatial resolution that can be 
achieved with diffusion-weighted spectroscopic measurements 
is similar to that of MRS methods in general, i.e., typical size 
of a volume of interest (VOI) in human experiments is several 
milliliters. 

Diffusion-weighted MRS (DWS) studies have been used to 
probe excised neural tissue (Assaf and Cohen, 1998a,b); In vivo 
DWS studies in small animals and primates were performed 
on muscle (Moonen et al., 1990; van Doom et al., 1996; de 
Graaf et al., 2000) and brain (Merboldt et al, 1993; van der 
Toorn et al, 1996; Pfeuffer et al, 2000; Dreher et al, 2001; 
Valette et al, 2007; Marchadour et al., 2012). Robust measure- 
ments have been performed in humans, (Kroenke et al, 2004; 
Ellegood et al, 2005a,b; Upadhyay et al, 2007, 2008). DWS 
studies at 7T provided tissue-specific information on several 
metabolites (Kan et al, 2012) and provided sensitive diffusion 
tensor MRS (DTS) data on NAA from small portion of the 
corpus callosum in a group of subjects with MS (Wood et al., 
2012). 

In this work we explore the potential of DWS to pro- 
vide compartment-specific information by focusing on a set of 
DWS results from a small portion of the corpus callosum. We 
show that the DWS data per se dependably reflects microstruc- 
tural features of the structures populated by the metabolites 
on which the measurements are performed. We then discuss 
the potential role that DWS may have in the evaluation of 
microstructural properties of tissue in conjunction with other 
methods, and its optimal positioning in the arsenal of MR-based 
microstructural tools. 

MATERIALS AND METHODS 
HUMAN SUBJECTS 

A total of 12 healthy volunteers (6 females, 6 males, age: 26 ± 4 
years) participated in this study. The study adhered to the LUMC 
Institutional Review Board guidelines and informed consent was 
obtained from all subjects prior to the study. 

MRI SCANNER/HARDWARE 

All experiments were performed on a 7 Tesla Philips Achieva 
whole-body MRI scanner (Philips Healthcare, Best, The 
Netherlands) equipped with gradient coils capable of a maxi- 
mum gradient strength of 40 mT/m and a slew rate of 200 T/m/s. 
A head coil consisting of a quadrature birdcage transmit and 
32-channel phased array receive was used for all measurements 
(Nova Medical Inc., Wilmington, MA, USA). 



MRI/MRS PROTOCOLS 
Anatomical images 

A short survey scan and a sensitivity encoding (SENSE) reference 
scan were followed by a 3D Ti -weighted gradient echo acquisition 
to allow for accurate planning of the experiment. Imaging 
parameters were: field of view (AP, FH, RL): 246.4 x 246.4 x 
174.0 mm 3 , resolution 1 x 1 x 1 mm 3 , TR/TE = 4.7/2.1 ms. Total 
scan time — 3 min. 

DTI protocol 

A single subject (male, 25 years) was scanned. Multi-slice single- 
shot 2D spin-echo echo planar imaging (EPI) (60 slices of 2 mm 
thickness, TR/TE = 9100ms/58 ms), field of view of 224 x 224 x 
120 mm 3 with an in-plane data matrix of 1 12 x 112, resulting in a 
data set with resolution of 2 x 2 x 2 mm 3 , 32 diffusion-weighting 
directions with b = 1000 s/mm 2 . Fat suppression was performed 
using spectral inversion recovery (SPIR). Parallel imaging was 
performed using SENSE with a reduction factor of 3 along the 
phase encoding direction (AP). Total scan time — 6 min. 

DWS 

The VOI for the experiment was the anterior body of the cor- 
pus callosum (Witelson segments A3-A4 and part of A2, see 
Figure 1 for a typical VOI positioning). VOI dimensions were 
25 mm(a-p), 10mm(r-l), and 8mm(f-h) for a total volume of 
2 cc. Slight variations in positioning of the VOI across subjects 
were necessary in order to adjust to slightly different patterns 
of anterior-posterior curvature. The sequence used for the DWS 
experiment was a point-resolved spectroscopy (PRESS) sequence 
(TE =121 ms) supplemented with a bipolar diffusion-weighting 
scheme (Kan et al., 2012). Cardiac synchronization on every 3rd 
cardiac cycle was achieved via a pulse-oximeter resulting in a TR 
of about 3 s and was used to avoid strong fluctuations in sig- 
nal intensity due to cardiac pulsation. Number of time-domain 
points = 2048, spectral width = 3000 Hz. Shimming up to second 
order was performed with pencil beam excitation and a typical 
tNAA line width was about 10 Hz throughout the experiment. 

Two directions with respect to the VOI coordinates were cho- 
sen for the diffusion weighting: (1) a pure right-left direction, 
thus mostly parallel to the direction of the callosal fibers); (2) one 
mostly perpendicular to the callosal fibers, forming a 45° angle 
between the anterior-posterior axis and the inferior-superior axis 
of the VOI (see diagram in Figure 1). These gradient directions 
can be denoted in the VOI coordinates as [1,0,0] and [0,1,1], and 




FIGURE 1 | Typical placement of the DWS VOI (A) on the sagittal plane 

and (B) on the coronal plane. The gradient directions used in the 
experiment are g[0,1 ,1 ], which is mostly perpendicular to the callosal fibers, 
and g[1 ,0,0], which is mostly parallel to them. 
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g[o,i,i] = \/2-g[i,o,o] for a specific choice of gradient strength. 
Diffusion-weighting parameters were: gradient pulse duration: 
34 ms, diffusion time (A): 60.5 ms, 8 equally spaced gradient 
amplitude values between 0 and 3.15 gauss/cm resulting in a set 
of fr-values of 0-6594 s/mm 2 for g[o,i,i] and 0-3297 s/mm 2 for 
£[i,o,o]- Each of the 15 direction/gradient amplitude conditions 
(2 directions x 7 gradient values + one acquisition without 
diffusion weighting) was acquired 48 times for a total of set of 
720 individual free induction decays (FIDs) and a total scan 
time of ~35 min. "Weak" water suppression was performed in 
order to preserve enough of the water signal (about 5-10 times 
the NAA singlet for all conditions) for subsequent phasing and 
frequency drift correction for each individual acquisition. To 
help enable this condition for the higher fo-values, water sup- 
pression was performed in a suboptimal manner by providing 
only a fraction of the RF power needed for full water suppres- 
sion. Following the DWS experiment, a short DWS acquisition 
(~3 min) with the same diffusion-weighting conditions and four 
acquisitions/condition was performed without water suppression 
and subsequently used for eddy current correction. 

Given the relatively large chemical shift displacement caused 
by the relatively low bandwidth of the 180° pulses in the PRESS 
sequence (about 4 mm for the chemical shift difference between 
the tNAA singlet (2.0 ppm) and the choline singlet (3.2 ppm) 
in the FH and LR directions), six of the subjects (4 females, 2 
males, 25.7 ± 4.8 years) were scanned with the frequency cen- 
tered around the tNAA peak and the other six (2 females, 4 males, 
27.3 ± 3.4 years) centered around the choline peak. For the six 
subjects scanned with the frequency at 2.0 ppm we report only 
the tNAA results, and for those scanned with the frequency at 
3.2 ppm, results are reported for both tCho and tCr, as the dis- 
placement of about 0.5 mm between the choline and the Cr/PCr 
singlet was deemed negligible. 



DWS DATA PROCESSING 

Data processing was performed using a custom written program 
in Matlab® (Mathworks Inc, Natick, MA, USA) as described pre- 
viously (Kan et al., 2012). The main steps in the data processing 
included: weighted summation and phasing of the individual out- 
puts of the 32 receive coils based on reference water signal, eddy 
current correction, zero-order phase correction, frequency drift 
correction, and subsequent averaging. The resulting spectra were 
analyzed using LCModel (Provencher, 1993) with an appropriate 
simulated set of basis spectra. 

MODELING OF THE DWS DIFFUSION DATA 

Analysis of DWS data has been performed in two ways. First, 
a non-biased analysis was carried out, where general diffusion 
properties for the three main metabolites were calculated: D\\ , D±, 
ADC, and pseudo-FA (D\\/D±). These properties were calculated 
based on the assumption of monoexponential decay with respect 
to fo-value: 



S(bi) = S(0) • e- b <-°< 



(1) 



where S(b) is the measured signal, S(0) is the signal without dif- 
fusion weighting, fo, is the fo-value in the direction i and D is the 
calculated diffusion coefficient in that direction. 

MODELED ANALYSIS OF THE DWS SIGNAL 

Expressions for signal attenuation of the NMR signal in the 
presence of magnetic field gradients, when the spin motion is 
restricted was first formulated by Neuman (1974) and applied 
to the Stejskal-Tanner experiment by Murday and Cotts (1968). 
The theory was extended to cover a variety of geometries (Balinov 
et al., 1994; van Gelderen et al., 1994) and succinctly formu- 
lated in the following expression (Aslund and Topgaard, 2009) for 
planar, cylindrical, and spherical pores: 



POST-PROCESSING AND ANALYSIS 

DTI processing was performed using ExploreDTI (Leemans et al, 
2009). A color-coded FA map and a 2D main eigenvector tick map 
were generated and displayed in Figure 2. 




R 
A 



FIGURE 2 | Color-coded FA map of a coronal slice through the body of 
the corpus callosum. The DWS VOI dimensions (20 mm RL x 8 mm FH) in 

that plane are represented by the yellow rectangle. The ticks in each voxel 
represent the direction of the main diffusion eigenvector. The resulting 
directional range is significant, approximately ±20°. 
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Where no = 2 for cylinders and 3 for spheres. a m satisfies the 
equation j[ (a m R) = 0, where j[ is the derivative of the Bessel 
function (for cylinders) or spherical Bessel function (for spheres) 
of the first kind and first order, and R is the radius of the cylinder 
or sphere, respectively. 

MODELING THE tNAA DWS SIGNAL 

Analysis of the tNAA signal was based on the assumption that 
in white matter tNAA resides exclusively in axons. An estimate 
of the typical axonal diameter in the callosal region in which the 
DWS measurements were performed was obtained from (Aboitiz 
et al., 1992). Accordingly, about 75% of axons in the region of 
interest possess an inner diameter such that D ax < 1.5 |xm, 
and for 24% of these axons 1.5 |xm < < 4.5 [im. These val- 
ues take into account about 30% post-mortem shrinkage in the 
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reported values, as described in Aboitiz et al. (1992). Axonal 
geometry inside the VOI is assumed to be a stack of curved cylin- 
ders, possessing a macroscopic curvature as well as a microscopic 
misalignment, or "fanning." In DTI analysis, macroscopic cur- 
vature of axonal tracts can be for the most part neglected, but 
in the case of DWS the dimensions of the VOI are of the same 
order of magnitude of the radius of curvature of the tract and thus 
the curvature cannot be neglected. This can be seen in Figure 2, 
in which a coronal slice of a directionally color-coded FA map 
is shown, along with a representation of a typical VOI in that 
plane. The variability of the fiber directionality within the VOI, as 
represented by the projections of the main diffusion tensor eigen- 
vectors on the plane, is clearly visible. Significant microscopic 
misalignment has been reported in the corpus callosum of the rat 
(Leergaard et al., 2010) and the macaque monkey (Sotiropoulos 
et al., 2012), and is reported in humans in this research topic by 
Budde et al., and thus has to be introduced as well. 

We base our formulation of the diffusion of tNAA in the 
axons on previous works (Assaf et al., 2004; Avram et al., 2008). 
For a single axon at an angle \|/ to the direction of the gradient, 
the diffusion-weighted signal S is the product of the diffusion 
perpendicular and parallel to the axonal wall: 

S (% A) = S ± (t ;± , A) ■ S,| {% , A) (3) 

where A is the diffusion time, ~qj = (8 is the duration of 
the pulsed field gradient, y is the gyromagnetic ratio of protons, 
g j is the diffusion- weighting gradient vector in the direction j), 

angle between and an axon in direction i. S± is given by 
Equation (1) expressed for the cylindrical case, and S is assumed 
to be governed by non-restricted diffusion, i.e.: 

S||(t*||, A) = S||(0)-e- t »- D (4) 

where b \\ = y 2 S 2 g 2 (A — |) and D is the diffusion coefficient 
of tNAA in the cytoplasm, from here onward D(tNAA). 

To account for both the macroscopic curvature within the 
VOI and the microscopic axonal misalignment, a distribution 
of the angle tyij is introduced, where = Qy + q>ij, where By- 
is the macroscopic contribution and (p y accounts for axonal 
misalignment. It is assumed that the displacement of NAA 
molecules during the diffusion time A is much smaller than 
the radius of curvature, hence 9,)- along this distance remains 
constant. In this case expression (3) is modified to yield the 
following form (Avram et al., 2008): 

N 

E P (^y) • sin (tyij) ■ S -L (<?; • c °s {^ij) , A) • 
;= i 

a \ S J (U ■ sin (jg) ; A ) 

S(<?;.A) = jj 

X>(*y) - sin (xj/y) 
l 

(5) 

For simplicity, it is assumed that the angles belong to a 
Gaussian distribution with zero mean and standard deviation a^, : 



1 =4- 

p (*|r|0, = 7=e 2 °* (6) 

ct^v 2jt 

A more appropriate description of angular distribution for 
the three-dimensional fanning of axons about a specific direc- 
tion is the Bingham distribution (Bingham and Mardia, 1978), 
applied in a number of studies that examined axonal fan- 
ning based on diffusion-weighted images (Kaden et al., 2007; 
Sotiropoulos et al., 2012). Here we chose the Gaussian distri- 
bution for simplicity, keeping in mind that the actual angular 
distribution, and in particular the portion attributed to the 
macroscopic curvature of the axonal tract, will be best described 
empirically rather than analytically. As the last step, the diffusion- 
weighted data in both directions was simultaneously fitted to the 
model in Equation (5) using two free-fitting variables: D(tNAA), 
the diffusion coefficient of tNAA in the direction parallel to 
the fibers, and 0^, the standard deviation of the total axonal 
dispersion. A diagram of the analysis procedures is given in 
Figure 3. 

MODELING THE tCr AND tCho DWS SIGNAL 

Analysis of the tCr and tCho signals must assume a significant 
contribution from glial content of these metabolites. The geo- 
metric representation of glia is rather complex, as the cells consist 
of a relatively large, oblong cell body with long, thin projections 
(Kuehnel, 2003; Verkhratsky and Butt, 2007). One extreme 
simplification is to assume that most of the cellular volume is 
represented by the cell body, approximated by a spherical pore. 
In order to obtain a representative compartment size, it is then 
possible to use Equation (2) for the spherical case. Since neither 
tCho nor tCr are exclusively astroglial, a complete analysis for 
Sm, the diffusion-weighted signal for the metabolite in question, 
should include an axonal fraction and a glial fraction, Pm,%\ and 
Pjvf.ax) respectively, weighing the DWS signal contribution from 
these compartments, Sjvr, gl and Sm,^'- 

S M [S, A, ~fd = P M , g i • S M , g i (r g i, D M , [§, A, ~gi]) 

+ ?M,ax • SM.ax^ax, D M , CT*, [8, A, g j]) (7) 

Here, Dm is the diffusion coefficient of the metabolite, r g i is the 
typical radius of the glial cell body (assumed to be a sphere), 
is the typical axonal diameter and cr^, is the axonal angular dis- 
persion. If cr^, is previously estimated through the DWS of the 
axonal constituent (tNAA), it can be used in fitting Equation 7 to 
the DWS data of the tCr or tCho, for further estimation of Pm,^, 
Pm, a*, and D m . 

STATISTICAL ANALYSIS OF RESIDUALS FOR THE MODELED DATA 

The quality of the fit of the experimental DWS data to the pro- 
posed models was assessed performing a j 2 goodness-of-fit test 
for the residuals. For this purpose, we used the Cramer-Rao 
Lower bounds (CRLB) of the spectroscopic measurements for the 
value of the variance of the single measurement. Typical CRLB 
values for the condition b = 0 s/mm 2 (in percent of the esti- 
mated peak integral) were 4% for the tNAA, 6% for the tCr, 
and 6% for the tCho. The highest CRLB values were those for 
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FIGURE 3 | A schematic representation of the modeling steps. In step 
(A) the tNAA DWS data in both diffusion-weighting directions is fitted to a 
model of diffusion in a set of misaligned cylinders with a gaussian angular 
distribution with standard deviation 04,. The fitting variables are the diffusion 
coefficient of tNAA and a^ . In step (B), the DWS signal of the two non-axonal 



specific metabolites, tCho and tCr is fitted to a mixed model of cylinders 
(axons) and spheres (glia), with cr^, obtained earlier as one of the input 
parameters. The fitting variables for each metabolite data are its diffusion 
coefficient, D(tCho) or D(tCr), and its glial fraction, P g i(tCho) or P g i(tCr), 
respectively. 



the maximum b-value in the [1,0,0] condition and typical val- 
ues were 18% for the tNAA, 13% for the tCr, and 13% for 
the tCho. 

RESULTS 

Typical diffusion-weighted spectra are shown in Figure 4. The 
spectrum in the top panel was acquired without diffusion weight- 
ing. The one below was acquired with the largest fo-value in the 
direction [0,1,1], which is mostly perpendicular to the callosal 
fibers. The bottom spectrum was acquired with b = 2422 s/mm 2 
in the [1,0,0] direction, which is roughly parallel to the callosal 
fibers. It is seen that the difference between the tNAA peak heights 
in the two diffusion-weighted spectra is significantly larger than 
the difference between the peak heights of the tCho and tCr. 
DWS data for the three metabolites: tNAA, tCho, and tCr, aver- 
aged over the six subjects scanned for tNAA and the six subjects 
scanned for tCho + tCr are shown in Figure 5. A qualitative esti- 
mate of the diffusion directionality, or anisotropy, for each of the 
three metabolites is provided by calculating D[i,o,o]/-D[0,i,i]> the 
ratio between a single diffusion coefficient calculated for the set 
of measurements of diffusion parallel, and one calculated for the 
diffusion perpendicular, to the fibers. Table 1 summarizes these 
data. It can be seen that D[i,o,o]/f [0,1,1] is the highest for tNAA 
(above 5), roughly twice the value for tCho and tCr. Interestingly, 
the values for Dro.ij] are rather similar for all three metabolites, 
D[o, 1,1] for tNAA being slightly lower than the other two, and sub- 
stantial differences are found between D[io,o] of tNAA and those 
of tCr and tCho. 



MODELING OF tNAA SIGNAL 

Figure 6 shows the DWS tNAA results from a single subject, to 
which the model described in the method section was fitted. The 
model fits the data well, well within the goodness-of-fit accep- 
tance criterion [x^j < 6.3, p > 0.9 for all subjects]. The two 
fitting variables were a^, the standard deviation of the axonal 
angular dispersion, and D(tNAA), the diffusion coefficient of 
tNAA in the cytoplasm. In order to exemplify the effect of ay, 
a suboptimal fit is shown in the left panel, in which cr^, is about 
50% of its optimal value. In that case, the diffusion effect is under- 
estimated for ^[0,1,1] (mostly perpendicular to the fibers) and 
overestimated for D[\ o,o] (mostly parallel to the fibers). The aver- 
age results for D(tNAA) and 0^ for the six subjects are given in 
Table 1. 

MODELING OF tCho AND tCr SIGNAL 

Choline soluble compounds are found in astrocytes and neu- 
ron/axons at a ratio of about 3:1, whereas Cr and PCr are 
ubiquitous in neural tissue, regardless of cell type (Choi et al, 
2007). A comprehensive analysis of the diffusion-weighted tCho 
and tCr signals should thus consider the contribution of both glial 
and axonal compartments, as is explicitly specified in Equation 7. 
In order to fit the data to the model in Equation 7, the axonal 
angular dispersion found for the tNAA was used as an input 
parameter for the modeling of the two remaining metabolite 
data, and the only fitting parameters were the glial fractions for 
the two metabolites, _P g i(tCho) and P g i(tCr), and the cytoplas- 
mic diffusion coefficients of tCho and tCr, D(tCho), and D(tCr), 
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FIGURE 4 | Typical spectra acquired from one of the subjects. (A) 

Spectrum acquired without diffusion weighting; (B) spectrum acquired in 
the direction [0,1,1] with b= 6594 s/mm 2 , the largest b-value used here; 
(C) spectrum acquired with b= 2422 s/mm 2 in the direction [1,0,0]. The 
strong diffusional anisotropy of the tNAA signal is clearly seen. 



respectively. The cell body radius was not a set as a free vari- 
able, but was gradually incremented until a stable fit with a 
constant error was obtained. This occurred at r g i = 9 (im, which 
is well within the range of astrocyte cell body radius (Kuehnel, 
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FIGURE 5 | Group-averaged data for the three main metabolites: 
top— tNAA; middle— tCr; bottom-tCho. 



2003). The model fits the data reasonably well, with x.( 14) < 11.7, 
p > 0.5 for all subjects. The group averages of the diffusion coeffi- 
cient and the glial fractions for tCho and tCr are given in Table 1 . 
The difference between P g i(tCr) and P g i(tCho) was statistically 
significant (p = 0.011, paired f-test). 

DISCUSSION 

The estimated D\\ and D± for the three metabolites measured 
in this work qualitatively display the most salient difference 
in compartmentation between tNAA, which is predominantly 
present in neuronal and axonal cytoplasm, and tCho and tCr, 
which are not unique to neurons and are also found in glial 
cells, mostly in astrocytes. In this work, the DWS data for tCr 
and tCho were acquired separately from those for tNAA, intro- 
ducing inter-individual differences as a possible confound, and 
it is thus desirable to devise a shortened acquisition proto- 
col that allows to acquire these data from the same subject. 
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Table 1 | Group-averaged DWS data (n = 6 for tNAA and n — 6 for tCr and tCho) for the three metabolites. 

ADC(mm 2 /s) D ± (mm 2 ) D (mm 2 /s) Dy/Dx D(model) (mm 2 /s) <r^,(°) P g | ia 

tNAA 0.24±0.01 0.079±0.003 0.41 ±0.01 5.2 ± 0.3 0.48±0.05 21.1 ±2.6 

tCr 0.18±0.01 0.095±0.004 0.26±0.02 2.8±0.4 0.37±0.06 - 0.4±0.1* 

tCho 0.16±0.02 0.093±0.007 0.22±0.02 2.4±0.4 0.42±0.08 - 0.5±0.1* 

*p = 0.077, paired t-test. 
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FIGURE 6 | tNAA DWS data from a single subject (symbols) fitted to 
the model described in the methods section and in step (A) in Figure 3. 

(A) The data fitted to the model with one fitting variable, D(tNAA), and 



= 10.0°; (B) the data fitted to the model with two fitting variables, 
D(tNAA), and cr^, . In this case D(tNAA) = 0.43 x 10~ 3 mm 2 /s and 
rj,j, = 22.3°. 



The good reproducibility of the DWS results for each metabo- 
lite, as can be appreciated from the results shown in Figure 5, 
is encouraging, and indicates that inter-subject variability is 
relatively low. From the structural perspective, it is useful to 
consider the intra-axonal space as cylinders with an inner diam- 
eter that is mostly in the 1-5 |xm range (Aboitiz et al, 1992), 
and astrocytes as represented by slightly oblong oval structures 
with average diameter of about 10-20 |xm (Glenn et al, 1992; 
Kuehnel, 2003) with thin processes that radiate radially from 
the soma in a rather isotropic fashion (Verkhratsky and Butt, 
2007). It has been established that tNAA is almost exclusively 
intracellular (Huang et al., 2000), and thus any contribution 
from "hindered" diffusion in the extracellular matrix can be 
neglected. The picture may be somewhat more complex for tCr 
and tCho, which both cross the blood brain barrier and are 
then cleared from the extracellular space through specialized 
transporters, although in normal conditions the levels of both 
extracellular tCr and tCho are low compared to intracellular lev- 
els (Klein et al., 1991; Loffelholz et al, 1993; Ohtsuki et al., 
2002). 

The strong structural anisotropy of axons is reflected by the 
significantly higher D[i o,o] value for tNAA than that for tCr 
and tCho. A strong incentive for an accurate estimation for 
D(tNAA) is the suggested link between axonal integrity and 



the diffusion coefficient along the axonal axis, or axial dif- 
fusivity as evaluated from the diffusion tensor (Song et al., 
2003; Kim et al., 2006; Budde et al, 2007). The axial diffu- 
sivity is the first eigenvalue of the diffusion tensor, Xi, and 
represents the diffusion coefficient in the main direction of 
the axonal tract. In DTI experiments, where the spatial res- 
olution is on the order of a few millimeters, Xi is mostly 
defined by two major factors: the "effective" diffusion coeffi- 
cient of water molecules in the intra- and extracellular space 
at the given diffusion time, and by the axonal organization, 
"fanning," or microscopic misalignment. The former is tightly 
related to the tortuosity of both the cytoplasmic as well as 
the interstitial medium, and it is hypothesized that structural 
damage to intra-axonal structures such as neurofilaments and 
microtubules resulting from pathological processes may affect 
the intra-axonal component of Xi. Accurate measurement of 
D(tNAA) has the potential to provide an additional source of 
information on intra-axonal structural integrity in white matter 
tracts, while avoiding the confound of the extra-axonal con- 
tribution. The potential sensitivity of D(tNAA) to pathology 
has been shown in the case of MS, where the axial diffusiv- 
ity of tNAA in the corpus callosum was significantly lower 
in a group of MS patients and correlated well with behav- 
ioral outcome (Wood et al, 2012). The specificity of DWS of 
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tNAA is offset by the low sensitivity of the method, which dic- 
tates significantly larger volumes than the typical DTI voxel and 
longer measurement times. One of the main problems resulting 
from using a large VOI is the impact of macroscopic curvature of 
the axonal tracts on the estimation of D (tNAA). In this work we 
showed that in the relatively simple case of the corpus callosum, it 
is possible to incorporate in the model a synthetic angular distri- 
bution of axonal directions. A simple and robust way to provide 
a realistic estimate for the macroscopic curvature can be based on 
an independent DTI data set, where the set of main eigenvectors 
for the DTI voxels within the DWS volume can be used to gener- 
ate a realistic distribution of angles between the tract directions 
and the diffusion-weighting gradients. These angular distribu- 
tions are patently non-Gaussian, and are expected to be different 
for each diffusion-weighting direction. We do not present such 
distributions in this work as we have not acquired DTI data 
from the subjects on which the reported DWS measurements 
were performed, but preliminary estimates based on data col- 
lected so far for a current project yield values of o$ [1,0,0] ~ 15° 
and 0e[O,l,l] ~ 10°, where 6 is the angle between the diffusion- 
weighting direction ([1,0,0] or [0,1,1]) and the main eigenvector 
of a DTI voxel within the DWS VOI. The quantitative incorpo- 
ration of these angular distributions in the processing pipeline of 
DWS data of NAA is one of the future directions we propose to 
investigate. Additionally, the addition of a DTI experiment to the 
DWS protocol will provide tractographic data for the visualiza- 
tion of white matter tracts of interest that are not as recognizable 
as the corpus callosum, and help position the DWS VOI on 
such tracts. 

The similar -Dpo.o] and D[o,ij] values for tCho and tCr ini- 
tially suggest that the compartmentation of these two metabolites 
in white matter is rather similar, and the differences can at least 
partly be explained by differences in molecular weight and in 
interactions with the surrounding medium. Upon more detailed 
analysis of the diffusion data for these two metabolites, a sig- 
nificant difference in their glial fraction is found, namely, the 
glial fraction of tCho for each subject is significantly higher than 
the glial fraction of tCr. This finding is supported by the liter- 
ature (Choi et al., 2007) and provides an encouraging support 
for the validity of the suggested model. It should be noted that 
since the measurements shown here were performed using a single 
diffusion time, it is impossible to use the data to extract infor- 
mation about compartment size, and in the proposed analysis 
the compartment size was either provided from histological data, 
as in the case of the axonal diameter, or given a lower limit in 
which the data fit well with the model, as in the glial diameter. 
The estimate of the latter (r g i = 9 |xm) agrees well with histo- 
logical data, but there is no doubt that measurements with a set 
of diffusion times would narrow down the estimates for D(tCr) 
and D(tCho), as well as those for the glial fraction for those 
metabolites. Moreover, we assumed here that the prevalent geom- 
etry of astrocytes is that of spherical compartments. This is a 
simplistic assumption that does not take into account astrocytic 
processes, which comprise a significant fraction of the astrocytic 
volume (Perge et al, 2009). The astrocytic processes are isotropi- 
cally distributed and thus do not affect the directional anisotropy 



of the DWS measurements, but since diffusion in fibers signif- 
icantly differs from diffusion in spheres in its dependence on 
diffusion time (Marchadour et al., 2012) it is expected that intro- 
ducing astrocytic processes to the model will have an impact 
on both the estimate of glial size as well as on D(tCr) and 
D(tCho). 

Diffusion properties of intracellular metabolites are not only 
indicative of cellular microstructure. In a recent study (Branzoli 
et al., 2013), DWS was applied to measure changes in ADC of 
metabolites in the primary visual cortex upon visual stimula- 
tion. In this work, the most salient result was that the ADC 
of tCr robustly increased with activation, and correlated well 
with the stimulation paradigm. Although this finding does not 
exclude the possibility of microscopic neuronal structural changes 
upon activation such as cell swelling (Le Bihan et al., 2006), 
the fact that the ADC of tNAA, a purely neuronal constituent, 
was significantly less correlated with the stimulus indicate a 
possible involvement of a metabolic process as responsible to 
the increase in ADC of tCr, rather than a structural one. The 
increase in SNR at ultra-high field facilitated the robust sep- 
arate assessment of ADC of glutamate and glutamine in gray 
and white matter (Kan et al, 2012). It is hoped that measur- 
ing the diffusion properties of glutamate and glutamine will 
provide further useful insights into such processes as glutamate 
release and the transport of glutamine from astrocytes to neurons 
upon neuronal activation, and contribute to the understanding 
of the role of glutamate compartmentation in disorders such as 
schizophrenia (Moghaddam et al., 1997; Moghaddam and Javitt, 
2012). 

In this work we provided some insight into the vast poten- 
tial of obtaining unique compartment-specific information from 
DWS, and an example of a scaffold for the analysis of such data. 
There are significant challenges associated with robust acquisi- 
tion, processing and analysis of diffusion-weighted spectroscopy. 
Additionally, the lower sensitivity of DWS when compared to 
diffusion-based imaging methods, prevents it from being used 
for high-resolution mapping of microstructural features of tis- 
sue. Conversely, the unique intrinsic compartment-specificity of 
the metabolites probed with DWS gives DWS a potential role 
in tissue characterization that water-based diffusion methods 
cannot fulfill. The challenge ahead is how to combine, in a 
multi-modal manner, DWS data with more sensitive but less- 
specific high-resolution data from e.g., DTI, to provide stronger 
explanatory power to both. Early attempts at doing so aimed 
at better understanding the processes that affect tissue in stroke 
(Wick et al., 1995; van der Toorn et al., 1996), and recently 
such attempt was aimed at characterizing tissue damage in nor- 
mal appearing white matter in MS (Wood et al, 2012). DWS 
will certainly gain from the advent of ultra-high field MRI for 
humans (7T and higher), as the sensitivity of localized spec- 
troscopy at higher fields is increased, with very few detrimental 
effects. It is hoped that the increase in availability of ultra-high 
field MRI scanners will generate an increase in interest in this 
unique tool, and that DWS will find its place in the growing arse- 
nal of methods used for the non-invasive investigation of tissue 
microstructure. 
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